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ABSTRACT 

We have undertaken a systematic search for candidate supernovae from high-redshift Population III 
stars in a field that has been observed with repeated imaging on a cadence of 2-3 weeks over a 2.2 
year baseline, the Spitzer /IRAC Dark Field. The individual epochs reach a typical 5 a depth of 1 ^Jy 
in IRAC Channel 1 (3.6 /mi). Requiring a minimum of four epochs coverage, the total effective area 
searched is 214arcminutes 2 . The unprecedented depth and multi-epochal nature of these data make 
it ideal for a first foray to detect transient objects which may be candidate luminous Pair Instability 
Supernovae from the primordial-metallicity first stars. The search was conducted over a broad range 
of timescales, allowing for different durations of the putative candidates' light curve plateau phases. 
All candidates were vetted by inspection of the Spitzer imaging data, as well as deep HST / ACS 
F814W imaging available over the full field. While many resolved-source objects were identified with 
Spitzer variability, no transient objects were identified that could plausibly be identified as high- 
redshift supernovae candidates. The resulting 95% confidence level upper limit is 23deg _2 yr _1 , for 
sources with plateau timescales under 400/(l+z) days and brightnesses above ~ 1 jtzJy. 

Subject headings: first stars:supernovae - population III stars: high redshift - infrared: IRAC 



1. INTRODUCTION 

Primordial metallicity Population III stars are thought 
to be the first luminous objects to form in the Uni- 
verse. Their formation marks the end of the cosmic 
dark ages and an important transition of the universe 
from a homogeneous state to a highly structured one. 
The UV photons produced by such stars at high red- 
shifts are also thought to be at least partly responsi- 
ble for re-ionizing the univers e (iTumlinson & Shul|2003, 
iBromm et alll200lL ISchaerenl2002L ISchaererll2003l ). It is 
believed that the explosive events that mark the end state 
of such st ars seed the interga l actic medium with heavy 
elements dGnedin fc Ostrikerl 119971 iFurlanetto fc Loea 
120031 iGreif et al.ll2007f K Hence, studying these objects 
is of great importance in helping us to understand the 
high-redshift universe. 

To date, the supernovae marking the deaths of first 
stars (Pop III SNe) at high redshift hay e not been ob- 
served (though see IWooslev et all [20071. We turn to 
theoretical modeling to gain a better understanding of 
the properties of such stars and when they might have 
existed. It is thought that Population III stars formed 
out of primordial-abundance H/He gas in low-mass dark 
matter halos. For primordial- abundance stars it is ex- 
pected that the explosion mechanism may drive not only 
"classical" supernovae, but also "hypernovae" for certain 
progenitor masses, driven through an electron-positron 
instability mechanism that results in explosive events 
with up to approximately one hund red times greater lu- 
minosities (jUmeda fe Nomotd 12002]) . Exactly when the 
epoch of the first stars began is still a ma tter of debate 
but e stimates place it at 10 < z < 50 (|Wise fe Abell 
2005). Thereafter, such objects could exist in primordial- 
metallicity pockets even at relatively low redshifts, even 



z < 2.5 (e.g. IScannapiecoll2005L iTornatore et al.ll2007h . 
The primordial metallicity of these stars leads to ineffi- 
cient cooling mechanisms through H2, leading to very 
high stellar masses and to a top-heavy initial mass 
function (IMF) with a large fraction of stars having 



M+ > lOOM^ (IBromm et all 
lAbel et all [2000 



19991 IBromm et alJl20u1 
20021 ) an d possibly even 



Abel et al.l . 

great er masses ( Bromm fe Loebll2004l lOmukai fc Pallal 
I2003T1 . 

In this paper we describe a search for candidate su- 
pernovae from high-redshift Population III stars in the 
Spitzer IRAC dark field. The target field is described 
in Section [2] In Section [3] we summarize the theoreti- 
cal modeling used as a guide for detecting transient ob- 
jects, which would include candidate Pop III SNe. Sec- 
tion [4] outlines the search methods used. In Section [5] 
we compare the Pop III SNe rate upper limits we find 
with expectations from the literature, and highlight how 
this type of search may be extended in the future. Where 
needed, we adopt H —73 kms -1 Mpc -1 , fl m — 0.27, and 
n A = 0.73. 

2. THE IRAC DARK FIELD 

The IRAC dark field is the dark current calib ration tar- 
get fo r the Infra-Red Array Camera (IRAC) (jFazio et al.l 
2004) on board the Spitzer space telescope. It is an ex- 
tragalactic field of very low background, in Spitzer's con- 
tinuous viewing zone near the North Ecliptic Pole. This 
area is observed at the start and end of each IRAC ob- 
serving campaign (2-3 weeks apart) since Spitzer's first 
light. For technical reasons anchored on the need for 
high-quality dark frames and on the normal precession 
of the observatory, there is only a modest overlap in the 
observed area on the sky from epoch to epoch, making 
for a position-dependent observed cadence in the time 
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domain data. 

The data used in this analysis are based on 128 dis- 
tinct epochs over the first 2.2 years of Spitzer's opera- 
tions. Each epoch is composed of multiple individual ex- 
posures at all of the available IRAC exposure times. The 
total full IRAC mosaic in the adopted dataset has a total 
observing time of > 500 hours, and covers an approxi- 
mately circular area, 20arcminutes in diameter. Each 
point in this mosaic typically has more than 10 hours of 
total integration time, with a maximum of ~ 80 hours in 
the area of maximal overlap across epochs. It is classi- 
cally confusion-limited in all four IRAC channels (3.6/im, 
4.5/!to, 5.8/im, S.Ofim). The I RAC source extract i on wa s 
done using Source Extractor, iBertin fc Arnoutsl (|1996l ) . 

The unprecedented depth and multi-epochal nature of 
these data make it ideal for a first foray into trying to 
detect supernovae from the first stars. The field will, in 
fact, continue to be observed through the full extent of 
the Spitzer Warm Mission, eventually giving a baseline 
of around seven (and possibly ten) years. There is now 
a wealth of multiwavelength data available in the dark 
field including Ch andra and HST / ACS F814W imaging 
(|Krick et al. l l2008l) . 

To find transient objects, a common method is to 
search for significant residuals when differencing regis- 
tered time-series observations. However, because the 
IRAC point spread function is asymmetric and fairly 
complex, and indeed has a different orientation on the 
image of each epoch, it is difficult to distinguish candi- 
date transient objects from artifacts in the difference im- 
ages. Therefore, the search is conducted through cross- 
correlating catalogs of objects detected in each individ- 
ual epoch. Each epoch of observation has a typical 5 a 
depth of 1/xJy, or m 3 . 6Aim (AB) ~ 23.9. The catalogs 
from all epochs were cross-matched with each other us- 
ing a 1 arcsec radius, resulting in a master database with 
31,492 sources, each of which had at least one detection 
in a single distinct epoch. This master database contains 
all the critical information for a transience search. For 
each object, the following are recorded: in which epochs 
the object is within the observed field and detected, and 
its flux density; in which epochs the object is within the 
observed field and not detected, and the known sensitiv- 
ity at that position (based on the integration time). 

If the spectral energy distributions of putative high- 
rcdshift Pop III SNe are above the measured detection 
limits for some period of time contained within the times- 
pan of this survey, they will appear as transient objects. 
We explore the expectations from theory and the poten- 
tial efficacy of our search in the next section. 

3. EXPECTATIONS FROM THEORY 

Massive Pop III stars with primordial metallicity are 
thought to be common at high redshift. Stars with 
M* < 140M© or M* > 26OM are thou ght to form 
black holes at the end of their evolution (e.g. iFrver et al.l 
120011 iHeeer fc Wooslevll2002D . Those which have masses 
between 140-260 M are thought to end their lives as 
pair- instability supernova (PISNe). Once helium burn- 
ing in the core of such stars has ceased the re is sufficient 
entro py to create positron-electron pairs ([Wise fc Abell 
2005). This process converts thermal energy to the mass 
of the particle pair and the pressure in the core is re- 
duced. This leads to a partial collapse which triggers 



TABLE 1 

Predicted Population III supernova rates found in the 

LITERATURE. 



dN/dz M progcn i tor z Reference 

(deg-y- 1 ) (a%) 

~0.2 250 ~ 10 Hcgcr ct ab (20021 ° 

0.34 100-500 10 Wise fc Abel (2005) 

50 250 > 15 Mackev et al. (2003) 

11 100 >13 Cen (2003)' 

25 140-260 5 Wcinmann fc Lilly (2005) 



a Rate of IHeeer e t_ab | ([2002|) incorpo rating correction factor de- 
termined by Wcinmann & Lilly (2005) and recalculated at z = 10 
'Corrected rate determined by Wcinmann & Lilly (2005) 



a thermonuclear explosion. The star is completely de- 
stroyed leading to a PISN, in which no remnant is left 
behind. At least one relatively local analog may al- 
ready have been observed (SN 2006 gy in NGC 1260; 
ISmith et al.ll2007l . IWooslev et al.ll2007l ). which lends sup- 
port to the possibility of this mechanism. PISNe would 
be "host-less" and as much as one hundred times more 
luminous than more-typical supernovae. 

To estimate the chances of detecting such objects, we 
turn to predictions from the literature for anticipated lu- 
minosities, durations (through thei r light curves), and 
frequency of events. Scannapicco (200 a), us ing light 
curves calculated by IWeinmann fc Lilly! (|2005l ) , predict 
peak apparent magnitudes of toab ~ 26.8 at z = 10 for 
250 Mq PISN, assuming negligible extinction. This sug- 
gests that the typical by-epoch depth of m 3 . 6Mm (AB) ~ 
24 of our search may be able to detect such objects at 
z~3-5. 

PIS N light curv e s are calc ulated in | Wise fc Abell 
(|2005h . iHeger et al.l (|2002f ). and IScannapiecol (|2005l ).A 
broad plateau phase is expected, which could last from 
~10 days to as long as a full year in the frame of the 
event. Since in the observed frame the light curve is 
stretched by a factor of (1+z), there could be events that 
would appear as near-continuous sources over the entire 
2.2-year monitoring span of the current dataset. As we 
ultimately anticipate a 7- to 10-year dataset by the end 
of the Spitzer Warm Mission, there is great potential 
in future analyses to encompass more plateau-duration 
possibilities. 

In Table Q] we list several predicted Pop III SNe differ- 
ential ra tes from the l iterat ure. The rates q uoted in Ta- 
ble ED for Eigiriiin (|2002l ) a ndlCenl (120031) incorporate 
the corrections determined by IWeinmann fc Lilly! (|2005l ) 
of (1 + z)~ 2 and (1 + z)' 1 , respectively. These rates 
are over redshift ranges largely beyond the sensitivity of 
our search, b ut are quoted here f or co mpleteness. By 
the review in Wcinma nn fc Lilly (2005), realistic rates 
are expected to be dN/d z ~ 4deg~ 2 yr~ 1 for z > 15 and 
0.2 deg-V^" 1 for z > 25. IWise fc Abel (EH) find a Pop 
III SNe rate of 0.34deg _2 yr _1 at z = 10, which changes 
negligibly over the mass range 1OOM < M* < 500M©. 
A wide range of values are expected, ~0.2— 50deg -2 
yr _1 , which is indicative of how the parameters involved 
in such prediction s are still not well-co nstrained. We also 
refer the reader to IScannapiecol (2005) for further discus- 
sion of these predictions. 
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4. SEARCH METHOD 

The search for transient objects is based on the mas- 
ter database described in Section [2] Systematic searches 
were conducted for three different ranges of possible 
plateau durations (in the observed frame), which are ap- 
propriate to the 2.2-year baseline of the present dataset, 
and which are plausible based on some of the theoret- 
ical expectations for the light curves. These are < 
t obs < 100, 100 < t obs < 200, and 200 < t obs < 400 
days. In addition to the practical convenience of ana- 
lyzing these data in this comparamentalized fashion, it 
also anticipates the possibility that some transient-source 
light curves may be truncated by either end of our 2.2 
year baseline by different timescales. 

We remind the reader that the light curves of all 
sources are largely discontiguous because the precise 
pointing (and orientation) varies epoch by epoch. This 
also affects the precise depth of each epoch, and changes 
the on-sky orientation of the IRAC diffraction pattern 
essentially leading to rotating diffraction spikes across 
images. We find that the photometry of objects near 
bright stars is systematically contaminated by these ro- 
tating diffraction spikes. Therefore we identify the stars 
and the sources close enough to them to be affected, 
and explicitly remove them from our master catalog. 
Imaging with HST '/ACS F814W was obtained between 
November & December 2006, one year after the end of 
the IRAC dataset described here. With this data, a 
plot of the Source Extractor (jBertin fe Arnoutsl Il996h 
ISOAREA value versus aperture magnitude photometry 
clearly separates the ridge of unresolved sources from re- 
solved galaxies, because unresolved sources have smaller 
ISOAREA than galaxies at any given magnitude. This 
method identifies 1147 unresolved sources. We empir- 
ically determine a radius of r = 0.06 x S^.Q^m arcsec 
(but no greater than 30 arcsec), to mask out a circular 
area around each star also removing an additional 5858 
sources which were in close proximity. 

Given the relative non-homogeneity in the survey's 
depth and cadence sampling, we apply several criteria 
to remove erroneous sources. 1) We impose an upper 
flux density limit of 40 /i Jy, which is an empirically de- 
termined practical threshold for removing additional ob- 
jects that are affected by diffraction spike artifacts, be- 
yond the masked radii, and object de-blending issues. 2) 
For all sources we only include epochs which have a 5 a 
detection of >l/xJy. 3) We do not consider sources that 
appear in only a single epoch. 4) Sources detected in 
only two or three epochs are only included in our search 
if they have 5 a detections of >1 /iJy in all epochs. 5) For 
any object with multiple epoch detections, we require the 
3.6fim photometric uncertainty to be less than 10% in at 
least half of its detected epochs. This ensures that this 
source is not rejected on the basis of large uncertainties 
in a small number of epochs (e.g. due to some epochs 
being particularly shallow relative to the others). 6) Fi- 
nally, we require at least two significant non-detections, 
to ensure a clear transient signal. 

After application of each of these criteria, 650 can- 
didates remained that warranted more careful follow-up. 
The 3.6/im light curves and the corresponding IRAC and 
HST imaging were all visually inspected by custom-built 
software that clearly identified all salient aspects for the 



object in the master database. No a priori restrictions 
were placed on the shape of the light curve. Several ob- 
jects with otherwise flat light curves were found to have 
a dramatic "flare-up" in a single epoch. Careful inspec- 
tion of the individual IRAC 100 s exposures of that epoch 
showed the spike to be a cosmic ray incident. The re- 
maining candidates were found to have resolved-source 
counterparts in the HST data, and are therefore can- 
didate low-redshift active galactic nuclei, which though 
extremely interesting in their own right, are clearly not 
Pop III SNe candidates. At the end of this careful analy- 
sis and vetting procedure, no viable candidates survived. 
In the next Section we calculate the formal limits implied 
by our search. 

5. DISCUSSION & CONCLUSIONS 

No Pop III SNe candidates to our sensitivity limit of 
7713.6^111 (AB) ~ 24 were identified. For the areal search 
over a total of 214arcminutes 2 , the rate limit is be- 
low 8deg~ 2 yr _1 . Using the same area but also mod- 
eling this non-detection limit as a Poisson distribution, 
the 95 %-confidence upper limit is 23deg~ 2 yr~ 1 , ignor- 
ing cosmic variance uncertainties. At z <[3, 5, 10] 
this also corresponds to a volumetric rate of [910, 480, 
270] Gpc" 3 yr" 1 , respectively. These limits are approx- 
imate, as the precise survey area relevant to each indi- 
vidual source is a complicated function of the varying 
field and depth by epoch. Furthermore, these limits only 
apply to moderate-duration events, with plateau phases 
lasting less than ~ 400/(1 + z ) days, by construction of 
the search. Many of th e predicted light curves given e.g. 
by iScannapiecol (|2005h have plateau phases up to 1 yr 
in the rest frame, which would be too long in the ob- 
served frame to be detected by the present survey, but 
which may be approachable in future incarnations of the 
Spitzer Dark Field su r veys. 

The iWise fc AbeH ([2005h and iHeger et alJ (pOOl ) 
z = 10 predicted differential rates of dN/dz ~ 
0.34deg~ 2 yr -1 and ^0.2deg~ 2 yr" 1 , respectively, are 
not ruled out, even if they appear at lower redshifts 
which may be above our sensitivity limit. Given that our 
search is mo st effectively probing z ~ 3 — 5, the differ- 
ential rate of lWeinmann fc Lilly! (|2005l ) of 25 deg -2 yr" 1 
at z = 5 is broadly comparable to our limit. It should 
be noted the very luminous PISNe that we would be sen- 
sitive to in this survey may be only a small fraction of 
all high-redshift supernova events. There is a distinc- 
tion currently being made between Pop III. 1 and III. 2 
stars, where the former class are of fully primordial abun- 
dance, and form in dark matter mini-halos, resulting in 
the pre-requisit e stellar masses of above ~ 100 to 
produ ce PISNe (| Johnson fc Brommll2006L iMcKee fc Tanl 
2008). These are distinct from the Pop III. 2 stars, which 
are expected to form through atomic c ooling processes, 
produ cing only ~ 10 M Q progenitors (jGreif fc Bromml 
l2006fh The PISNe Pop III.l pro genitors may consist o f 
only some 10% of Pop III SNe (jGreif fc Bromml [20061) . 
Furthermore, the "pristine" Pop III.l p rogenitors can 
suffer dramatic negative feedback (e.g. IMcKee fc Tanl 
2008), which may additionally limit their relative num- 
bers. These we re all consider a tions not yet taken in the 
predictions by iMackev et al.l ((2003), which led to the 
high expectation rates. 

To do a comprehensive search for Pop III SNe using 
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Spitzer (or a future platform with similar capabilities) 
would require a survey area of 1 deg 2 with exposures of 
5000 seconds (50 x 100 s exposures) per point on the sky. 
This would reach ttiab ~ 26 at the 3-5 a level, depending 
on extinction. Each epoch would require approximately 
200 hours of observation to cover the survey area with a 5 
arcminute field of view. This greater depth complicates 
the IRAC data reduction as confusion is a significant, 
but tractable problem. Trading depth for greater area 
may not be optimal for detecting the /iig/iesf-redshift Pop 
III SNc, which are expected to be quite faint, ttiab <^ 
26. This requirement is increasingly relaxed for lower 
redshift events, which may possibly exist at redshifts as 
low as z ~ 2. This ideal survey would need to be imaged 
every two months for a period of several years in order 
to plausibly sample much or all of the plateau phase of 
the light curve. The total program would therefore be of 
the order of 2000-3000 hours. 

In the meanwhile, the current dataset continues to 
expand as Spitzer continues to operate, acquiring new 
IRAC observations every 2-3 weeks. Including obser- 
vations made during the Spitzer Warm Mission, the full 
dataset will span at least seven years, and maybe as many 
as ten full years. Future searches following the procedure 
described here should be powerful for identifying plausi- 
ble candidates, or at least in setting firmer limits on the 
production rate of PISNe at high redshift, and thus set- 
ting practical limits on the relative abundances of Pop 
III. 1 versus Pop III. 2 progenitors ultimately informing 
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